Comment on 'Effect of on-site Coulomb repulsion on superconductivity in the 
boson-fermion model' (Phys. Rev. B66, 134512 (2002)) 
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The two-dimensional (2D) boson-fermion model (BFM) of high-temperature superconductors, 
numerically studied by Domanski (Phys. Rev. B66, 134512 (2002)) is not a superconductor. The 
critical temperature of the model is zero for any symmetry of the order parameter. The opposite 
conclusion advocated by Domanski stems from a mean-field approximation (MFA) neglecting the 
boson self-energy which is qualitatively erroneous for any-dimensional BFM. 

PACS numbers: PACS: 71.20.-z,74.20.Mn, 74.20.Rp, 74.25. Dw 



Recently Domaski p| and some other authors (see ref- 
erences in pj) claimed that 2D BFM with hybridised 
fermions and immobile hard-core bosons is capable to 
reproduce the phase diagram of cuprates. The model is 
defined by the Hamiltonian, 

H = E ^L^ + EoJ2 b %+ (!) 

k,<T=T,l q 

gN~ 1/2 ^2 (0k& q C-k+q/2,T C k+q/2,i + H - C ) I 
q,k 

where £k = — 2t(cosfca: + cosfcj,) — a is the 2D energy 
spectrum of fermions, Eq = A B — 2u is the bare bo- 
son energy with respect to the chemical potential //, g 
is the magnitude of the anisotropic hybridisation inter- 
action, 0k = (/>_k is the anisotropy factor, and N is the 
number of cells. Ref. argued that 'superconductivity 
is induced in this model from the anisotropic charge ex- 
change interaction (g<^k) between the conduction-band 
fermions and the immobile hard-core bosons', and 'the 
on-site Coulomb repulsion U competes with this pair- 
ing' reducing the critical temperature T c less than by 
25%. The author of Ref. Q neglected our study of BFM 
, which revealed a devastating effect of the boson self- 
energy on T c . Here I show that because of this effect 
T c = OK in the model, Eq.(l), even in the absence of the 
Coulomb repulsion, U = 0, and the mean-field approx- 
imation of Ref. [l| is meaningless for any-dimensional 
BFM. 

Using MFA Ref. decouples bosons and fermions in 
Eq.(l) replacing boson operators by c-numbers for q = 0. 
Then T c is numerically calculated using the linearised 
BCS equation for the fermionic order-parameter Ak, 



E N 



A k > tanh(&'/2fc B T) 



(2) 



where Eq is determined by the atomic density of bosons 
{n B ) as (Eq.(9) in Ref. Q) 
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En 
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(3) 



with T = T c and g 2 = g 2 (l - 2n B ). While Eq.(2) is 
perfectly correct, Eq.(3) is incorrect because the boson 



self-energy Sf,(q, O n ) due to the exchange interaction is 
missing. At first sight Q the self-energy is small if g 
is small in comparison to the kinetic energy of fermions. 
However Eb(0, 0) diverges logarithmically at zero temper- 
ature 2], no matter how week the interaction is. There- 
fore it should be kept in the density sum-rule, Eq.(3). 
Introducing the boson Green's function 



Diet, n r , 



1 - 2n L 



£ -£fc(q, fir. 



(4) 



one must replace incorrect Eq.(3) by 



k B T 
N 



(5) 



where r = +0, and n n = 2irk B Tn (n = 0,±1,±2...). 
The divergent (cooperon) contribution to S&(q, O n ) is 
given by 



£&(qA) = -|^E^ x 



(6) 



tanh[g k _ q/a /(2A: B r)] + tanhE k+q/a /(2A: B r)] 

Ck-q/2 + £k+q/2 — 

so that one obtains 



E b (q,0) =S fe (0,0) 



2M* 



0(q 4 



(7) 



for small q with any anisotropy factor compatible with 
the point-group symmetry of the cuprates. Here M* is 
the boson mass, calculated analytically in Ref.0] with 
the isotropic exchange interaction and parabolic fermion 
band dispersion (see also Ref. 3]), and h = 1. The BCS- 
like equation (2) has a nontrivial solution for Ak, if 



E a = -E 6 (0,0). 



(8) 



Substituting Eq.(7) and Eq.(8) into the sum-rule, Eq.(5) 
one obtains a logarithmically divergent integral with re- 
spect to q, and 



const =Q _ 



Jo d ll<l 



(9) 
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Using a 'bubble' approximation for the self-energy Ref.0 
proved that the Cooper pairing of fcrmions in BFM is im- 
possible without the Bose-Einstein condensation (BEC) 
of real bosons. The bubble approximation is actually ex- 
act because of the logarithmic divergence of the Cooperon 
diagram, as was also confirmed in Ref.jj]. Hence, the 
devastating result, Eq.(9) is a direct consequence of the 
well-known theorem, which states that BEC is impossible 
in 2D. 

One may erroneously believe that MFA results pj can 
be still applied in three-dimensions, where BEC is pos- 
sible. However, increasing dimensionality does not make 
MFA a meaningful approximation for the boson-fermion 
model. This approximation leads to a naive conclusion 
that a BCS-like superconducting state occurs below the 
critical temperature T c ~ p exp (— Eq / z c ) via fermion 
pairs being virtually excited into unoccupied bosonic 
states H Iff. Here z c = g 2 N(0) and N(0) is the density 
of states (DOS) in the fermionic band near the Fermi 
level fi. However, the Cooper pairing of fermions is im- 
possible via virtual unoccupied bosonic states. It occurs 
only simultaneously with the Bose-Einstein condensation 
of real bosons in the exact theory of 3D BFM . 

The origin of the simultaneous condensation of the 
fermionic and bosonic fields in 3D BFM lies in the soften- 
ing of the boson mode at T = T c caused by its hybridiza- 
tion with fermions. Indeed, Eq.(8) does not depend on 
the dimensionality, so that the analytical continuation 
of Eq.(4) to real frequencies u> yields the partial boson 
DOS as p(uj) = (1 - 2n B )S(uj) at T = T c for q = in 
any-dimensional BFM for any coupling with fermions. 

Taking into account the boson damping and dispersion 
shows that the boson spectrum significantly changes for 
all momenta. Continuing the self-energy, Eq.(6) to real 
frequencies yields the damping (i.e. the imaginary part 
of the self-energy) as 



cosh(gf + w/(4fc B T c )) 



cosh(-g£ + u/{Ak B T c )) 



(10) 



where £ = vp/ (4fcsT c ) is a coherence length. The damp- 
ing is significant when q£ << 1. In this region 7(q, ui) — 
LUTTz c / (8fcsT c ) is comparable or even larger than the bo- 
son energy oj. Hence bosons look like overdamped dif- 
fusive modes, rather than quasiparticles in the long- 
wave limit 0, 1 contrary to the erroneous conclusion of 
Ref.0, that there is 'the onset of coherent free-particle- 
like motion of the bosons' in this limit. Only outside the 
long-wave region, the damping becomes small. Indeed, 
using Eq.(10) one obtains 7(q, <jS) — lwkz c / (2qvp) « u>, 
so that bosons at q » l/£ are well defined quasiparti- 
cles with a logarithmic dispersion, uo(q) = z c ln(g£) @. 
As a result the boson dispersion is distributed over the 
whole energy interval from zero up to Eq, but not a delta- 
function at Eq even in the weak-coupling limit. 

The main mathematical problem with MFA stems from 
the density sum rule, Eq.(5) which determines the chem- 
ical potential of the system and consequently the bare 
boson energy Eq(T) as a function of temperature. In 



the framework of MFA one takes the bare boson en- 
ergy in Eq.(2) as a temperature independent parameter, 
Eq = z c \n(fi/T c ), or determines it from the conserva- 
tion of the total number of particles, Eq.(5) neglecting 
the boson self-energy £&(q,fi) 0, H H). Then Eq.(2) 
looks like the conventional mean-field BCS equation, or 
the Ginzburg-Landau equation (near the transition) with 
a negative coefficient a oc T — T c at T < T c in the linear 
term with respect to A(T). Hence, one concludes that 
the phase transition is almost the conventional BCS-like 
transition, at least at Eq >■ T c [a, E|, and, using the 
Gor'kov expansion in powers of A, finds a finite upper 
critical field H C 2(T) These findings are mathemati- 
cally and physically wrong. Indeed, the term of the sum 
in Eq.(5) with Q„ = is given by the integral 



T 



f/q 



27r 3 S + S b (q,0) 



(11) 



The integral converges, if and only if Eq ^ — £{,(0, 0). In 
fact, Eq + Efe(Q, 0) is strictly zero in the Bose-condensed 
state, because \ib = — \Eq + £&(0, 0)] corresponds to the 
boson chemical potential relative to the lower edge of the 
boson energy spectrum. More generally, fj,b = corre- 
sponds to the appearance of the Goldstone-Bogoliubov 
mode due to a broken symmetry below T c . This exact 
result makes the BSC equation (2) simply an identity Q 
with a = at any temperature below T c . On the other 
hand, MFA violates the density sum-rule, predicting the 
wrong negative a(T) below T c . 

Since a(T) = 0, the Levanyuk-Ginzburg parameter || 
is infinite, Gi — oo. It means that the phase transi- 
tion is never a BCS-like second-order phase transition 
even at large Eq and small g. In fact, the transition is 
driven by the Bose-Einstein condensation of real bosons 
with q — 0, which occur due to the complete softening 
of their spectrum at T c in 3D BFM. Remarkably, the 
conventional upper critical field, determined as the field, 
where a non-trivial solution of the linearised Gor'kov 
equation occurs, is zero in BFM, H C 2(T) = 0, because 
a(T) = below T c . It is not a finite H c2 (T) found in 
Ref. Q using MFA. Even at temperatures well below 
T c the condensed state is fundamentally different from 
the BCS-like MFA ground state, because of the pairing 
of bosons 10]. The pair-boson condensate significantly 
modifies the thermodynamic properties of the condensed 
BFM compared with the MFA predictions. 

This qualitative failure of MFA might be rather un- 
expected, if one believes that bosons in Eq.(l) play the 
same role as phonons in the BCS superconductor. This 
is not the case for two reasons. The first one is the den- 
sity sum- rule, Eq.(5), for bosons which is not applied to 
phonons. The second being that the boson self-energy 
is given by the divergent (at T = 0) Cooperon diagram, 
while the self-energy of phonons is finite at small cou- 
pling. 

I have to conclude that the numerical work by Doman- 
ski [l| does not make any sense in any dimension. There 
is nothing in 2D BFM to compete with because the model 
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is not a superconductor even without the Coulomb re- 
pulsion. The MFA results || do not make any sense 
in three dimensions either, because the divergent self- 
energy has been neglected in calculating T c and H C 2(T). 



The common wisdom that at weak coupling the boson- 
fermion model is adequately described by the BCS the- 
ory, is negated by our results. 
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